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Forecasting  F10.7  with  solar  magnetic  flux 
transport  modeling 
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[i]  A  new  method  is  presented  here  to  forecast  the  solar  10.7  cm  (2.8  GHz)  radio  flux,  abbreviated  F10.7, 
utilizing  advanced  predictions  of  the  global  solar  magnetic  field  generated  by  a  flux  transport  model. 
Using  indices  derived  from  the  absolute  value  of  the  solar  magnetic  field,  we  find  good  correlation 
between  the  observed  photospheric  magnetic  activity  and  the  observed  Fi0.7  values.  Comparing 
magnetogram  data  observed  within  6  hours  of  the  F10  7  measurements  during  the  years  1993  through 
2010,  the  Spearman  correlation  coefficient,  rs,  for  an  empirical  model  of  Fi0  7  is  found  to  be  0.98.  In 
addition,  we  find  little  change  in  the  empirical  model  coefficients  and  correlations  between  the  first  and 
second  9  year  intervals  of  the  18  year  period  investigated.  By  evolving  solar  magnetic  synoptic  maps 
forward  1-7  days,  this  new  method  provides  a  realistic  estimation  of  the  Earth-side  solar  magnetic 
field  distribution  used  to  forecast  F10.7.  Spearman  correlation  values  of  approximately  0.97,  0.95,  and 
0.93  are  found  for  1  day,  3  day,  and  7  day  forecasts,  respectively.  The  method  presented  here  can  be 
expanded  to  forecast  other  space  weather  parameters,  e.g.,  total  solar  irradiance  and  extreme 
ultraviolet  flux.  In  addition,  near-term  improvements  to  the  F10.7  forecasting  method,  e.g.,  including 
far-side  magnetic  data  with  solar  magnetic  flux  transport,  are  discussed. 

Citation:  Henney,  C.  J.,  W.  A.  Toussaint,  S.  M.  White,  and  C.  N.  Arge  (2012),  Forecasting  F10  7  with  solar  magnetic 
flux  transport  modeling.  Space  Weather,  10,  S02011,  doi:10.1029/2011SW000748. 


1.  Introduction 

[2]  Forecasting  variations  in  the  solar  ultraviolet  (UV) 
and  extreme  ultraviolet  (EUV)  flux  is  of  great  interest  for 
space  weather.  To  a  large  extent  these  fluxes  control  con¬ 
ditions  in  the  ionosphere  that  play  a  major  role  in  com¬ 
munications  and  navigation  signals,  e.g.,  of  the  ionospheric 
total  electron  content  [e.g.,  Richards  et  ah,  1994;  Maruyama, 
2010].  Current  UV  and  EUV  forecast  models  commonly 
rely  on  auto-regressive  and  time  series  analysis  of  past 
solar  measurements  of,  e.g.,  the  sunspot  number  (SSN)  and 
solar  10.7  cm  (2.8  GHz)  radio  flux  [e.g.,  Cliatterjee,  2001;  Lean 
et  ah,  2009;  Huang  et  ah,  2009],  commonly  abbreviated  as 
F10.7.  Parameters  such  as  SSN  and  F10.7  are,  in  general, 
proxies  for  solar  magnetic  activity.  For  example,  the  SSN  is 
commonly  used  to  estimate  F10.7  for  a  given  day.  Though 
the  SSN  and  F10.7  relationship  has  correlated  well  in  the 
past,  Svalgaard  and  Hudson  [2010]  found  this  may  no  longer 
be  the  case  for  the  current  solar  cycle.  The  correlation 
change  with  cycle  may  be  linked  to  the  observation  that  the 
magnetic  field  strengths  of  sunspot  umbrae  appear  to  be 


1Space  Vehicles  Directorate,  Air  Force  Research  Laboratory, 
Kirtland  Air  Force  Base,  New  Mexico,  USA. 

2National  Solar  Observatory,  Tucson,  Arizona,  USA. 


varying  with  time  [Penn  and  Livingston,  2006;  Livingston  and 
Penn,  2009].  These  findings  highlight  the  need  to  look  for 
better  magnetic  proxies  for  solar  activity.  We  show  here 
that  the  observed  Earth-side  solar  magnetic  field  strength 
and  distribution  can  be  used  to  estimate  F10.7  surprisingly 
well. 

[3]  Chapman  and  Boyden  [1986]  demonstrated  that  pho¬ 
tospheric  magnetic  field  measurements  can  be  used  to 
model  the  solar  total  irradiance.  By  using  disk-position 
weighted  sums  of  weak  (10  <  IBI  <  100  G)  and  strong  (IBI  > 
100  G)  magnetic  fields  (B)  separately.  Chapman  and  Boy¬ 
den  estimated  the  fluctuations  of  facular  and  sunspot 
regions.  Using  the  same  magnetic  sum  limits,  without  the 
disk-position  weighting,  Ulrich  [1991]  and  Parker  et  al. 
[1998]  revealed  that  such  magnetic  field  indices  could  be 
used  as  a  proxy  for  F10  7.  These  sums  of  the  absolute  value 
of  the  solar  magnetic  field,  measured  from  the  150  foot 
Solar  Tower  at  the  Mount  Wilson  Observatory  (MWO),  are 
referred  to  as  the  MPSI  (Magnetic  Plage  Strength  Index) 
and  MWSI  (Mount  Wilson  Spot  Index).  Extending  the 
period  of  study  by  Chapman  and  Boyden  to  three  solar 
cycles,  Jain  and  Hasan  [2004]  also  found  that  the  MPSI  and 
MWSI  indices  can  be  used  to  model  the  solar  irradiance, 
while  Schmahl  and  Kundu  [1995,  1998]  showed  that  radio 
measurements  at  a  number  of  frequencies  can  be  used  to 
reproduce  solar  irradiance  variations. 
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[4]  Expanding  on  work  of  Ulrich  [1991]  and  Parker  et  al. 
[1998],  we  outline  here  a  method  for  forecasting  space 
weather-related  solar  indices  using  future  estimates  of 
the  global  flux  distribution  on  the  photosphere.  Utilizing 
global  photospheric  magnetic  maps  of  the  solar  surface, 
the  future  Earth-side  magnetic  activity  is  estimated  from  a 
flux  transport  model.  We  have  selected  F10  7  for  this  study 
since  it  is  an  extremely  accurate  and  reliable  data  product 
widely  used  as  a  proxy  for  solar  magnetic  activity.  The 
analysis  presented  here  could  also  be  expanded  to  forecast 
EUV/UV  and  total  solar  irradiance  [e.g.,  Wenzler  et  al, 
2006;  Krivova  et  al,  2009].  Solar  sources  of  Fi0.7  and  mag¬ 
netic  flux  transport  models,  along  with  the  method  used 
here  to  forecast  Fi0.7  and  preliminary  results,  are  discussed 
in  this  paper. 

2.  Solar  Sources  of  F10.7 

[5]  Solar  10.7  cm  radio  flux  data  have  become  one  of  the 
most  widely  used  space  weather  indices  for  a  number  of 
reasons:  F107  arises  at  least  in  part  from  the  same  hot 
coronal  plasma  that  produces  the  Sun's  ionizing  (UV,  EUV 
and  soft  X-ray)  flux  [e.g..  Swamp  et  al,  1963;  Tapping  and 
DeTracey,  1990];  the  measurement  of  F10.7  is  carried  out 
using  a  very  careful  and  highly  repeatable  calibration 
procedure  [Tapping  and  Cliarrois,  1994];  and  there  is  an 
important  historical  record  of  measurements  back  to  1947 
that  can  be  used,  and  is  used  here,  to  investigate  long-term 
trends  and  correlations  with  other  phenomena  [Covington 
andMedd,  1949;  Covington,  1969;  Tapping,  1987]. 

[6]  The  emitting  wavelength  of  10.7  cm  is  at  a  critical 
point  in  the  Sun's  radio  spectrum  [e.g.,  Schmahl  and  Kundu, 
1998;  White,  1999].  At  longer  wavelengths,  e.g.,  20  cm,  the 
variable  component  of  the  Sun's  radio  emission  is  domi¬ 
nated  by  thermal  bremsstrahlung  from  the  hot  gas  in  the 
solar  corona,  which  can  be  optically  thick.  At  shorter 
wavelengths,  e.g.,  6  cm,  the  corona  is  always  optically  thin 
to  thermal  bremsstrahlung  but  regions  of  strong  magnetic 
field  in  the  corona  can  be  optically  thick  due  to  gyroreso- 
nance  emission  [e.g..  White  and  Kundu,  1997].  Where  the 
ambient  magnetic  fields  are  strong  enough  for  the  electron 
gyrofrequency  to  be  a  low  harmonic  of  the  observing  fre¬ 
quency,  thermal  gyroresonance  becomes  an  important 
emission  process.  In  weaker  magnetic  fields,  thermal 
bremsstrahlung  emission  continues  to  be  dominant.  The 
wavelength  of  10.7  cm  lies  in  a  regime  where  we  expect 
thermal  bremsstrahlung  to  be  largely  optically  thin,  and 
thus  large  gyroresonance  sources  over  active  regions  can 
be  seen  through  the  coronal  plasma  as  bright  contributors 
to  the  radio  flux.  Plages  outside  active  regions  will  con¬ 
tribute  to  Fi0.7  via  thermal  bremsstrahlung,  but,  unlike 
active  regions,  plages  generally  do  not  have  field  strengths 
sufficiently  strong  in  the  corona  to  contribute  gyroreso¬ 
nance  emission  to  Fi0.7.  The  relative  contributions  of 
bremsstrahlung  and  gyroresonance  emission  to  F10.7 
remain  a  matter  of  contention:  Felli  et  al  [1981]  and 
Schmahl  and  Kundu  [1998]  argue  that  gyroresonance  emis¬ 
sion  dominates  the  rotationally  modulated  component  of 


F10.7,  while  Tapping  and  DeTracey  [1990]  argue  that  thermal 
bremsstrahlung  is  the  main  contributor.  Unfortunately,  to 
date  there  have  been  no  systematic  high-spatial-resolution 
imaging  measurements  at  the  wavelength  of  10.7  cm  that 
can  settle  this  question. 

[7]  Whichever  emission  dominates  the  variable  compo¬ 
nent  of  F10  7,  we  expect  to  observe  at  least  a  qualitative 
relationship  between  Fi0.7  and  magnetic  activity.  The 
emission  measure  of  material  at  coronal  temperatures 
over  an  active  region  that  produces  thermal  brems¬ 
strahlung  at  radio  wavelengths  and  line  emission  at  EUV 
wavelengths  is  believed  to  reflect  the  amount  of  coronal 
heating  taking  place  on  magnetic  field  lines:  heat  depos¬ 
ited  in  the  corona  is  thought  to  propagate  down  field  lines 
and  heat  chromospheric  material  at  the  foot  points  that  in 
turn  rises  back  into  the  corona  [e.g.,  Cargill  and  Klimchuk, 
2004,  and  references  therein].  In  turn,  the  rate  of  coronal 
heating  is  found  to  be  linked  to  the  amount  of  magnetic 
flux  in  the  active  region  [e.g.,  Golub  et  al,  1980].  Therefore 
the  thermal  bremsstrahlung  emission  at  radio  wave¬ 
lengths  depends  strongly  on  magnetic  features  at  the 
photosphere,  and  it  is  clear  that  gyroresonance  emission, 
which  requires  strong  magnetic  fields  in  the  solar  corona, 
must  depend  on  the  magnetic  field  in  active  regions.  While 
a  few  previous  studies  have  addressed  this  relationship 
[Ulrich,  1991;  Parker  et  al,  1998],  there  has  been  extensive 
work  comparing  F107  with  the  sunspot  number  [e.g., 
Bouwer,  1992;  Floyd  et  al,  2005;  Tapping  et  al,  2007;  Svalgaard 
and  Hudson,  2010;  Johnson,  2011;  Tapping  and  Valdes,  2011]. 
In  this  paper  we  determine  a  quantitative  relationship 
between  F10.7  and  solar  magnetic  fields  and  then  take 
advantage  of  state-of-the-art  modeling  of  future  magnetic 
field  distributions  to  predict  Fi0.7  several  days  ahead. 

3.  Global  Solar  Magnetic  Flux  Transport 

[8]  Solar  magnetic  flux  transport  models  provide  an 
estimate  of  the  instantaneous  global  spatial  distribution  of 
the  solar  magnetic  field.  The  goal  of  these  models  is  to 
evolve  the  synoptic  magnetic  flux  distribution  and  match 
observations  by  incorporating  rotational,  meridional,  and 
supergranular  diffusive  surface  transport  processes  to 
predict  the  magnetic  field  in  locations  where  direct  mea¬ 
surements  are  not  available  [e.g.,  Devore  et  al,  1984;  Sheeley 
et  al,  1987;  Worden  and  Harvey,  2000;  Schrijver  and  De  Rosa, 
2003,  and  references  therein]. 

[9]  One  of  the  primary  uses  of  such  instantaneous  global 
maps  is  as  input  to  coronal  and  solar  wind  forecasting 
models,  e.g.,  the  Wang-Sheeley-Arge  (WSA)  model  [Arge 
and  Pizzo,  2000;  Arge  et  al,  2003,  2004].  Output  from  mag¬ 
netic  flux  transport  models  aid  to  minimize  non-physical 
monopole  moments  in  the  magnetic  maps  used  to  initi¬ 
alize  the  coronal  and  solar  wind  models.  Such  monopoles 
commonly  occur  in  Carrington  synoptic  maps  near  the 
eastern  limb  as  a  result  of  merging  newly  observed  data 
and  during  periods  when  the  solar  polar  regions  are  not 
well  observed  from  the  Earth.  The  magnetic  flux  transport 
model  used  here  is  the  Air  Force  Data  Assimilation 
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Table  1.  Empirical  F10  7  Model  Coefficients  (m0,  in  u  and  m2  From  Equation  (1))  and  Correlation  Coefficients  From  the  Regression 
Fits  for  the  18  Year  Period  Studied  (1993-2010),  Along  With  the  Corresponding  Standard  Deviations2 


Period 

Year  Range 

Coverage 

<Forecast  Time> 
(days) 

m0  (sfu) 

m1  (sfu/G) 

m2  (sfu/G) 

rs 

r 

Total 

1993-2010 

67.0% 

0.13  ±0.05 

66.08  ±0.09 

8.508  ±0.11 

16.56  ±0.18 

0.98 

0.97 

1.13  ±0.05 

65.26  ±0.10 

8.522  ±0.11 

17.20  ±0.19 

0.97 

0.96 

3.13  ±0.05 

65.00  ±0.10 

10.76  ±0.10 

15.90  ±0.22 

0.96 

0.95 

7.13  ±0.05 

64.75  ±0.09 

14.83  ±0.09 

8.861  ±0.27 

0.94 

0.91 

A 

1993-2001 

70.9% 

0.15  ±0.05 

65.41  ±0.14 

8.153  ±0.14 

16.81  ±0.22 

0.98 

0.97 

1.15  ±0.05 

64.88  ±0.15 

8.446  ±0.14 

17.04  ±0.24 

0.98 

0.96 

3.15  ±0.05 

64.80  ±0.15 

10.90  ±0.13 

15.46  ±0.27 

0.96 

0.94 

7.15  ±0.05 

64.59  ±0.14 

15.56  ±0.13 

6.870  ±0.34 

0.94 

0.90 

B 

2002-2010 

63.2% 

0.12  ±0.05 

66.31  ±0.12 

9.058  ±0.18 

16.92  ±0.31 

0.97 

0.98 

1.12  ±0.05 

65.44  ±0.13 

8.495  ±0.17 

18.22  ±0.33 

0.96 

0.97 

3.12  ±0.05 

65.14  ±0.13 

10.50  ±0.16 

17.02  ±0.37 

0.94 

0.96 

7.12  ±0.05 

64.92  ±0.12 

13.87  ±0.14 

11.95  ±0.46 

0.92 

0.93 

“The  fit  and  correlation  values  are  provided  for  the  period  split  into  two  9  year  intervals:  A  (1993-2001)  and  B  (2002-2010).  The  correlation 
values  include  the  Spearman  (rs)  and  Pearson  (r)  coefficients.  The  percentage  of  days  used  within  the  period  ranges  (i.e..  Total,  A,  and  B)  are 
shown  in  the  coverage  column.  The  table  rows,  within  each  of  the  three  ranges,  are  for  the  mean  forecast  times.  Note  that  the  forecast  times 
are  slightly  greater  than  the  target  time  (i.e.,  0,  1,  3,  and  7  days). 


Photospheric  Flux  Transport  model,  referred  to  as  ADAPT  the  last  map  was  created.  The  choice  of  KPVT  and  VSM 

[Arge  et  al,  2010,  2011].  ADAPT  utilizes  a  modified  version  data  instead  of,  e.g.,  MDI  (Michelson  Doppler  Inter- 

of  the  Worden  and  Harvey  (WH)  flux  transport  model  ferometer)  or  MWO,  is  arbitrary  and  based  on  the  devel- 

[Wordeti  and  Harvey,  2000].  opment  version  of  the  ADAPT  code.  Other  magnetograph 

[10]  The  modified  WH  flux  transport  model  within  input  can  be  used  with  the  method  outlined  below,  with 

ADAPT  employs  an  ensemble  of  model  realizations  based  the  caveat  that  the  inferred  photospheric  field  strengths 

on  using  a  range  of  different  model  parameters  con-  between  instruments  can  vary  as  much  as  a  factor  of  2 

strained  by  the  estimated  errors  of  each  parameter  (e.g.,  [Jones  et  al,  1993;  Wenzler  et  al,  2004;  Thornton  and  Jones, 

meridional  drift).  The  data  assimilation  method  currently  2002;  Riley,  2007].  Jones  et  al.  [2004]  found  that  the  mag- 

incorporated  by  ADAPT  is  the  ensemble  least  squares  netic  fields  measured  by  the  KPVT  and  VSM  agree,  gen- 

method  [Arge  et  al,  2011].  The  use  of  an  ensemble  of  erally,  to  within  a  few  percent.  The  VSM  magnetic 

models  is  beneficial  in  two  significant  ways:  first,  the  sensitivity  improved  with  new  cameras  in  late  2009.  In 

ensemble  variance  is  utilized  to  estimate  the  model  error  total,  for  the  18  years  studied  here,  there  are  two  magnetic 

when  determining  the  relative  weight  given  to  new  data  strength  discontinuities  on  the  order  of  a  few  percent: 

and  model  calculations  in  the  assimilation  step;  second,  September  2003  due  to  the  switch  from  the  KPVT  to  the 

comparing  the  different  ensemble  members  individually  VSM  and  January  2010  due  to  a  camera  upgrade  in  the 

with  observations  (e.g.,  solar  polar  field  strength  and  cor-  VSM. 

onal  hole  boundaries)  provides  feedback  to  determine  the  [12]  For  this  study  the  ADAPT  code  was  modified  to 

magnitude  and  time  dependent  behavior  of  large  scale  create  forecast  maps  0  to  7  days  in  the  future  at  20  UT  to 

magnetic  flow  patterns  such  as  the  meridional  drift  rates  in  synchronize  with  the  F10  7  measurements  near  local  noon 

each  hemisphere.  The  F10  7  forecast  method  outlined  in  the  at  Penticton,  British  Columbia,  Canada  (earlier  measure¬ 
following  section  uses  forecast  synoptic  maps  generated  ments  were  made  at  Ottawa,  Ontario  [Tapping,  1987]).  The 

by  ADAPT  to  provide  estimates  of  the  Earth-side  solar  0  day  maps  are  created  to  model  how  well  nearly  simul- 
magnetic  field  distribution  1  to  7  days  into  the  future.  taneous  photospheric  magnetic  observations  correlate 

with  measured  values  of  F10  7.  The  0  day  map  is  generated 
for  20  UT  using  the  most  recent  available  observed  mag- 
4.  Forecasting  F10.7  With  Flux  Transport  netogram  (note  that  KPVT  and  VSM  data  are  not  acquired 

[11]  The  ADAPT  global  photospheric  maps  used  for  this  at  the  same  time  each  <%)•  Since  the  most  recent  data  are 

study  were  created  using  line-of-sight  magnetogram  data  always  acquired  before  20  UT,  and  ADAPT  then  advances 

from  1993  through  2010  from  the  Kitt  Peak  Vacuum  Tele-  the  assimilated  data  to  20  UT,  there  is  variability  in  the 

scope  (KPVT  [Jones  et  al,  1992])  and  Vector  Spectro-  forecast  time.  The  analysis  here  is  restricted  to  days  on 

magnetograph  (VSM  [Henney  et  al,  2009]).  A  new  map  is  which  new  data  are  available  within  a  6  hour  window 

generated  by  ADAPT  each  time  an  observed  magneto-  before  20  UT‘  With  this  acceptance  window,  the  0  day 

gram  is  available.  Typically,  the  KPVT  and  VSM  full-disk  maPs  used  for  this  study  have  been  evolved  forward 

magnetograms  are  available  at  a  cadence  of  approximately  3A  hours'  on  average,  and  thus  will  not  perfectly  represent 

one  per  day.  If  no  magnetograms  are  available,  a  new  photospheric  magnetic  conditions  at  20  UT  when  F10  7  is 

ADAPT  map  is  generated  at  a  cadence  of  ~24  hours  since  measured.  The  fourth  column  of  Table  1  provides  the 
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Figure  1.  The  weighted  sums  of  the  absolute  magnetic  fields,  from  Air  Force  Data  Assimila¬ 
tion  Photospheric  Flux  Transport  model  (ADAPT)  Earth-side  pixels,  with  values  greater  than 
25  gauss  and  less  than  150  gauss  (i.e.,  a  plage  index,  SP,  similar  to  the  Magnetic  Plage  Strength 
Index  (MPSI))  and  greater  than  or  equal  to  150  gauss  (i.e.,  an  active  region  index,  Sa,  similar  to 
the  Mount  Wilson  Spot  Index  (MWSI)). 


mean  and  standard  deviation  of  the  departure  of  the  KPVT 
and  VSM  magnetogram  observation  time  from  20  UT  for 
the  entire  18  years,  and  the  first  and  last  nine  year  periods. 

[13]  Comparable  to  the  empirical  model  used  by 
Chapman  and  Boyden  [1986]  and  Jain  and  Hasan  [2004]  for 
modeling  the  total  solar  irradiance,  we  estimate  the  F10.7 
flux  based  on  a  fit  to  ADAPT-generated  magnetic  fields 
using  ranges  designed  to  separate  the  plage-dominated 
component  from  the  sunspot-dominated  component, 
analogous  to  the  Magnetic  Plage  Strength  Index  and  the 
Mount  Wilson  Spot  Index.  We  tested  over  a  dozen  other 
magnetic  field  ranges  and  combinations,  but  found  the  fit 
and  correlation  improvements  to  be  negligible.  In  addi¬ 
tion,  subdividing  the  magnetic  fields  into  more  than  three 
ranges  resulted  in  alternating  negative  coefficients.  Since 
the  fit  improvements  were  negligible  and  do  not  justify  the 
increased  degrees  of  freedom,  we  chose  to  use  two  mag¬ 
netic  ranges. 

[14]  From  ADAPT  Earth-side  map  pixels,  we  compute 
the  sums  of  the  absolute  magnetic  field  values  greater  than 
25  gauss  and  less  than  150  gauss  (i.e.,  a  plage  index,  SP, 
similar  to  MPSI)  and  greater  than  or  equal  to  150  gauss 
(i.e.,  an  active  region  index,  Sa,  similar  to  MWSI),  to  model 
Fi0.7.  Since  the  synoptic  maps  produced  by  ADAPT  are  in 
heliographic  coordinates  and  hence  essentially  indepen¬ 
dent  of  the  Sun-Earth  distance,  we  model  the  values  of 
Fio.7  adjusted  to  1 AU  in  solar  flux  units  (sfu),  where  1  sfu  is 
equivalent  to  10  22  W  m  2  Hz-1.  The  empirical  model 
used  here  is  expressed  as: 

Fmodel  =  WO  +  OTl  SP  +  m2SA  ( 1 ) 

where 


^  IBrl<150G 

SP  =  = —  V  IB  rlu# 

UJe  25G<IBrl 

Sa  =  IB  rUe. 

Z^  150G<IBrl 

The  variable  ug  is  a  surface  area  weighting  relative  to 
latitude  to  convert  the  ADAPT  map  (i.e.,  180  latitude  by 
360  longitude  pixels)  values  to  equal  area.  The  normal¬ 
ization  sum,  identical  for  both  SP  and  Sa,  is  over  all  of  the 
1802  Earth-side  pixels.  The  magnetic  field  is  represented 
by  Br  to  emphasize  that  the  fields  are  radial  for  the  ADAPT 
heliographic  frame  maps,  where  the  radial  field  is  derived 
from  the  magnetograms  assuming  that  the  observed  field 
values  are  the  line-of-sight  component  of  a  radial  field.  We 
derived  the  best-fit  coefficients  for  equation  (1)  separately 
for  0,  1,  3,  and  7  day  ADAPT  forecasts  generated  from  all 
available  daily  magnetograms,  for  the  total  18  years  and 
the  first  and  last  nine  years,  which  are  shown  in  Table  1. 

[15]  Note  that  the  default  1802  pixel  area  used  encom¬ 
passes  pixels  that  would  not  normally  be  observed  in  a 
sky-frame  full-disk  image.  Earlier  versions  of  the  forecast 
method  estimated  the  actual  observed  area  as  viewed  from 
Earth,  taking  into  account  the  solar  tilt  angle  relative  to 
Earth's  orbit,  along  with  the  viewing  angle  relative  to  the 
Earth-Sun  distance.  The  differences  in  spatial  weighting, 
however,  resulted  in  negligible  improvements  in  the  final 
fits  and  correlations.  With  the  low  resolution  of  the 
ADAPT  maps,  and  the  lack  of  the  strong  fields  in  polar 
regions  that  are  normally  associated  with  F10.7,  we  chose  to 
use  1802  area  from  the  maps  for  this  preliminary  study. 
The  sum  values  for  both  Sp  and  Sa  are  shown  in  Figure  1. 
In  the  figure,  note  that  the  SA  values  are  rather  flat  during 
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Figure  2.  Adjusted  observed  F10.7  (solid  line)  and  the  empirical  model  (dotted  line), 
smoothed  with  a  21  day  running  mean,  for  the  18  year  period  studied.  Note  that  periods  with 
straight  lines  are  generally  associated  with  Vector  Spectromagnetograph  (VSM)  instrument 
downtime  (e.g.,  April-June  2004,  February-April  2006,  and  December  2009). 


the  exceptional  low  activity  period  during  late  2008  and 
early  2009  compared  to  the  previous  minimum.  The  higher 
variability  of  the  stronger  field  regions  is  most  likely  a 
result  of  the  smaller  spatial  scale  and  shorter  lifetime  of 
strong  fields,  i.e.,  sunspot  regions. 

5.  Results  and  Discussion 

[i6]  The  0  day  ADAPT-derived  F10.7  values  and  the 
observed  (adjusted  to  1  AU)  F10.7  values  for  the  full  18  year 
period  (1993  through  2010)  are  compared  in  Figure  2.  The 
figure  highlights  the  good  agreement  between  the  two  time 
series,  but  also  reveals  that  the  model  values  are  generally 
too  large  during  minimum  solar  magnetic  activity  and  too 
weak  during  solar  maximum  activity.  A  more  detailed 
comparison  between  the  observed  (forecast  valid-time) 


values  and  the  1  day  and  3  day  forecasts  is  exhibited  in 
Figure  3.  With  new  activity  rotating  onto  the  Earth-side  for 
the  first  time,  the  ADAPT  maps  are  expected  to  exhibit  a 
phase  lag  since  they  do  not  incorporate  solar  far-side  data, 
such  as  helioseismic  inferences  of  the  emergence  of  new 
regions.  Without  information  about  activity  on  the  far  side 
of  the  Sun,  the  first  knowledge  of  new  magnetic  activity  is 
its  appearance  on  the  east  limb  provided  by  solar  rotation, 
where  observational  errors  are  typically  much  higher  than 
the  model  errors.  The  observational  errors  arise  primarily 
from  foreshortening  and  the  canopy  effect,  along  with  the 
highly  variable  horizontal  magnetic  signal  that  increases 
toward  the  limb  [Harvey  et  ah,  2007].  Since  the  model  error 
is  small  by  comparison,  the  observed  data  can  take  a  few 
days  to  be  fully  assimilated. 


Date  (day-month,  year) 

Figure  3.  An  example  period  of  increasing  solar  magnetic  activity  showing  daily  values  of  the 
adjusted  observed  F10.7  (solid  line),  modeled  1  day  (diamond),  and  3  day  (plus  sign)  forecast 
estimations. 
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Figure  4.  Interval  A  (1993-2001),  the  percent  difference  between  the  0  day  (solid),  3  day 
(dotted),  and  7  day  (dashed)  ADAPT  forecasts,  and  the  adjusted  observed  F10.7  values. 
Smoothed  with  a  21  day  running  mean,  the  differences  highlight  how  the  forecasts  generally 
underestimate  F10.7  during  high  activity  (1999-2001). 


[17]  Figures  4  and  5  display  the  percent  difference 
between  the  0,  3,  and  7  day  model  values  and  observations 
for  the  18  years  split  into  two  intervals:  1993  through  2001 
(interval  A)  and  2002  through  2010  (interval  B).  These  fig¬ 
ures  highlight  when  the  empirical  model  is  best  (during 
low  activity  variability)  and  worst  (during  high  variability) 
with  differences  between  the  model  and  observations 
ranging  from  a  few  percent  to  25%  for  the  3  day  values. 
Besides  the  lack  of  far-side  magnetic  activity  information, 
this  may  reflect  the  fact  that  the  strong  field  regions  are  not 
well  represented  with  the  nominal  spatial  resolution  of 
ADAPT  maps,  i.e.,  360  by  180  (longitude  by  latitude)  pixels. 
It  may  also  represent  a  change  in  the  relative  contributions 
to  F10.7  by  gyroresonance  and  thermal  bremsstrahlung 
emission  at  high  activity  levels.  Comparing  Figures  1 
and  2,  note  that  nearly  all  the  variability  of  the  Fm7  signal 


correlates  well  with  the  magnetic  active  region  sum,  SA. 
Attempts  to  fit  the  F10.7  signal  with  only  SA  resulted  in  a 
good  fit  in  general,  however,  the  local  minima  were 
underestimated.  That  is,  once  prominent  activity  rotated 
off  the  Earth-side,  the  F10.7  background  level  is  not  well 
matched  with  SA  alone,  and  SP,  representing  weaker  fields, 
is  needed  to  encompass  the  observed  variability. 

[is]  For  the  3  day  ADAPT  forecast  maps,  nearly  80%  of 
the  Earth-side  pixels  have  been  evolved  for  72  hours  or 
longer  and  ~20%  of  the  pixels  have  been  evolved  for 
approximately  13  days  or  greater  since  the  last  Earth-side 
magnetogram  observation.  If  the  flux  transport  model  of 
ADAPT  were  "perfect,"  the  0  day  Fi0.7  empirical  model 
coefficients  would  be  expected  to  be  equally  valid  for  the  1 
through  7  day  forecast  maps,  however,  the  mid-range 
magnetic  fields  tend  to  decrease  faster  in  the  model 


Figure  5.  Interval  B  (2002-2010),  the  percent  difference  between  the  0  day  (solid),  3  day 
(dotted),  and  7  day  (dashed)  ADAPT  forecasts,  and  the  adjusted  observed  F107  values. 
Smoothed  with  a  21  day  running  mean,  the  differences  highlight  how  the  forecasts  generally 
underestimate  F10.7  during  high  activity  (2002-2006). 
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calculations  than  in  the  observations.  The  change  is 
reflected  in  the  increase  of  the  rti\  coefficient  value  from 
0  day  to  7  day  forecasts  for  all  period  fits  in  Table  1.  The 
general  amplitude  drop  could  be  the  result  of  new 
background  (i.e.,  weak  field)  and  active  region  emergence 
not  accounted  for  or  too  much  flux  diffusion  comprised 
within  the  modeled  supergranular  flow.  Since  the  Worden 
and  Harvey  flux  transport  within  the  ADAPT  model 
applies  far  less  diffusion  to  the  strong  magnetic  fields  [see 
Worden  and  Harvey,  2000],  the  drop  is  expected  to  be  pri¬ 
marily  from  the  weak  fields.  By  fitting  the  different  fore¬ 
cast  maps  (e.g.,  0  day  and  7  day)  separately,  this  lack  of 
mid-range  flux  can  clearly  be  accounted  for  by  the  quality 
of  the  correlations.  The  mid-range  flux  drop  was  not 
detected  until  this  work,  and  highlights  how  forecasting 
global  indices  can  provide  unique  feedback  to  flux  trans¬ 
port  models.  The  lack  of  additional  emergence  or  exces¬ 
sive  diffusion  will  be  addressed  by  an  ongoing  study  to 
validate  the  flux  transport  model  on  varying  timescales, 
from  12  hours  to  6  days  with  Earth-side  activity  and  13  to 
15  days  for  returning  activity. 

[19]  The  difference  between  the  fit  coefficients  for  the 
two  9  year  intervals  and  the  total  period  is  found  to  be 
small,  as  depicted  in  Table  1,  and  well  within  the  uncer¬ 
tainty  in  the  scaling  between  the  KPVT  and  VSM  magne¬ 
tograms  for  this  period.  The  effect  of  this  difference  is 
minimal  as  demonstrated  by  the  fact  that  using  the  inter¬ 
val  A  coefficients  to  determine  the  interval  B  radio  flux 
values  results  in  a  correlation  decrease  of  less  than  one 
half  of  a  percent.  The  Spearman  correlation  values,  using 
interval  A  coefficients,  match  those  found  fitting  interval  B 
independently  with  0.97,  0.96,  0.94,  and  0.94  for  the  0, 1,  3, 
and  7  day  forecasts  respectively.  The  ability  to  interchange 
coefficients  supports  the  argument  that  the  relationship 
between  F10.7  and  the  SA,p  parameters  does  not  exhibit 
cycle-to-cycle  variations.  In  addition,  the  stability  of  the 
coefficients  reinforces  the  idea  that  magnetic  parameters 
similar  to  SA,p  are  to  be  preferred  as  solar  activity  proxies 
over  parameters  such  as  the  sunspot  number,  which  has 
recently  been  shown  to  have  a  cycle-to-cycle  dependence 
[Penn  and  Livingston,  2006;  Livingston  and  Penn,  2009].  When 
the  18  year  period  is  subdivided  into  spans  less  than 
9  years  (e.g.,  5  years),  the  fits  are  dominated  by  either  high 
activity  or  inactivity  leading  to  poor  forecasting  once  the 
activity  trend  changes.  In  addition,  as  stated  earlier,  we 
chose  to  use  two  magnetic  ranges  for  the  model  even 
though  the  fits  could  be  improved  (albeit  very  slightly)  by 
splitting  the  magnetic  range  into  three  or  more  bins.  We 
found  that  the  coefficients  from  fits  using  more  than  two 
ranges  differed  greatly  between  the  two  intervals,  most 
likely  because  the  additional  bins  simply  allow  for  more 
degrees  of  freedom  for  the  fit. 

6.  Near-Future  Goals  and  Improvements 

[20]  A  better  understanding  of  the  sources  of  the  solar 
10.7  cm  radio  flux  signal  is  critical  to  improving  the  current 
empirical  model  development  and  evolving  to  a  model 


atmospheric  solution.  For  example,  higher  field  thresholds 
for  the  active  region  index  in  equation  (1)  were  tested  with 
negligible  improvement  and  greater  uncertainty  in  the 
fitted  coefficients.  With  higher  spatial  resolution  maps, 
however,  we  might  expect  improved  fits  using  a  higher 
field  threshold  corresponding  to  the  third  harmonic  of  the 
gyrofrequency  [e.g..  White  and  Kundu,  1997],  i.e.,  333  gauss 
in  the  corona,  once  the  relative  contribution  of  brems- 
strahlung  and  gyroresonance  emission  is  better  under¬ 
stood.  With  full-disk  F107  observations  from  the  EVLA 
(Expanded  Very  Large  Array)  proposed  in  the  near-future, 
we  expect  to  achieve  a  more  detailed  understanding  of  the 
F10  7  source  locations  in  relation  to  the  photospheric  mag¬ 
netic  field.  The  planned  comparison  work  between  EVLA 
and  magnetogram  image  data  is  anticipated  to  aid  in  the 
determination  of  the  minimum  spatial  resolution  for  the 
global  magnetic  maps  needed  to  improve  future  F10.7 
forecast  models.  In  addition,  the  observational  study  is 
expected  to  provide  feedback  on  the  limitations  of  a  simple 
empirical  model  outlined  in  the  previous  section.  The 
benefit  of  using  magnetic  flux  transport  with  global  solar 
magnetic  maps  is  that  we  can  move  away  from  an 
empirical  model  toward  a  physics-based  model  for  fore¬ 
casting  F10.7  [e.g.,  Lionello  et  al.,  2009],  Ultimately,  daily 
spatially  resolved  F10  7  observations  [see  Saint-Hilaire  et  al., 
2012]  would  permit  long-term  feedback  and  refinement  of 
flux  transport  models  like  ADAPT. 

[21]  While  the  current  version  of  the  ADAPT-driven 
empirical  model  proves  to  be  practical  for  estimating  the 
F10.7  value  for  forecast  timescales  of  1  to  3  days,  including 
solar  far-side  magnetic  activity  estimates  will  be  key  to 
modeling  more  accurate  forecast  values  and  extending 
forecasts  to  7  days.  The  far-side  estimates  are  possible  with 
the  helioseismic  acoustic  holography  technique  [e.g., 
Lindsey  and  Brann,  1990,  1997]  that  can  detect  large  mag¬ 
netic  regions  on  the  solar  far-side.  The  helioseismic  far- 
side  detection  technique  has  recently  been  parametrized 
in  terms  of  photospheric  magnetic  field  strength  [Gonzalez 
Hernandez  et  al.,  2007],  allowing  for  a  quantitative  estima¬ 
tion  of  the  solar  magnetic  activity  on  the  solar  far-side 
while  updating  the  Earth-side  of  global  synoptic  maps. 
The  practical  application  of  the  far-side  data  has  recently 
been  demonstrated  to  agree  well  with  the  observed  Fm7 
signal  [Gonzalez  Hernandez  etal.,  2011].  The  inclusion  of  far- 
side  data  with  ADAPT  is  part  of  a  newly  funded  3  year 
study.  The  ADAPT  model  is  ideal  for  representing  the 
nominal  distribution  and  morphology  of  the  far-side  active 
region  polarity  (e.g.,  Hale's  rule)  and  orientation  (e.g.,  Joy's 
rule)  within  the  ensemble  realizations.  Though  the  flux 
polarity  and  distribution  will  have  large  uncertainty,  a 
good  estimate  of  the  inferred  magnetic  flux  location  and 
magnitude  is  all  that  is  required  to  improve  the  F10.7 
forecast. 

7.  Summary 

[22]  The  solar  10.7  cm  radio  flux  forecast  method  pre¬ 
sented  here  is  shown  to  agree  well  with  observation. 


7  of  9 


Approved  for  public  release;  distribution  is  unlimited. 


S02011 


HENNEY  ET  AL.:  FORECASTING  F10  7  WITH  FLUX  TRANSPORT 


S02011 


resulting  in  Spearman  correlation  values  of  0.97,  0.95  and 
0.93  for  1,  3,  and  7  day  forecasts,  respectively.  The 
empirical  model  developed  here  is  based  on  the  estima¬ 
tion  of  the  future  global  solar  photospheric  magnetic  field 
using  the  flux  transport  of  the  ADAPT  model.  By  evolving 
solar  magnetic  synoptic  maps  forward  1  to  7  days,  this 
method  provides  a  realistic  estimation  of  the  Earth-side 
solar  magnetic  field  distribution  and  a  key  preliminary 
step  toward  more  physics-driven  forecasting  of  space 
weather  parameters.  As  the  ADAPT  model  is  updated  in 
the  near-future  with  solar  far-side  input,  along  with  flux 
transport  model  parameter  improvements  from  forth¬ 
coming  validation  work,  we  expect  improved  agreement 
between  forecasts  and  F10.7  observations.  Ultimately,  the 
disagreement  between  forecast  values  and  observations 
may  provide  useful  feedback  for  the  selection  of  ADAPT 
model  parameters  within  the  model  ensemble.  In  addi¬ 
tion,  with  higher  spatial  resolution  we  fully  expect  that  flux 
transport  models  similar  to  ADAPT  could  be  used  to 
forecast  the  total  solar  irradiance  and  the  sunspot  number, 
in  addition  to  solar  wind,  Fm7,  ultraviolet  and  extreme 
ultraviolet  flux. 
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